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Abstract 

Possible variations in the dynamical behaviour of L1YF4 due to its structural changes follow- 
ing several pressure-induced phase transitions are examined by making use of the complementary 
techniques of quasi-harmonic lattice dynamics and molecular dynamics simulation. The phonon 
spectra in the entire Brillouin zone and the respective Gibbs free energies are calculated for the 
three high-pressure polymorphs of L1YF4 that are stable at T = 0, with an aim to understand 
their relative stabilities as functions of pressure and temperature in terms of volume compression 
and vibrational entropy. Molecular dynamics simulations provide qualitative impressions about 
a temperature-driven second-order transformation and also of kinetic effects in the subsequent 

pressure-driven first-order phase transition. In addition, the calculations predict anomalous ther- 
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mal expansion at low temperature in phases I and Ha while irreversibilty of phase II — ► phase III 
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transition on subsequent pressure release. 



I. INTRODUCTION 

Ternary halide LiYF 4 (a laser host material of current interest |l|, |2j, |3, |4J , H |y, [D, |s|, |3, |l0|) 
displays several structural phase transitions J4J, |y, |8[ upon compression which turn out to be 
unique in comparison of isostructural CaW0 4 . For example, while the latter has a scheel- 
ite (J4i/o) to wolframite (P2\/c) kind of high-pressure phase transformation, the former 
shows an intermediate fergusonite phase (12 /a, above room temperature; P2i/c, in the low 
temperature region), associated with a soft phonon mode.|5| Recently, the post-fergusonite 
phase of LiYF 4 has been predicted as wolframitelike monoclinic one.J8|,|9( Interestingly, the 
axial ratio in both the compounds evolves in a different way under the effect of pressure. |l0|| 
In view of the above, we go for a comparative analysis of various dynamical properties that 
these phases (Fig. ^) are likely to demonstrate. To reach this, we carry out extensive lat- 
tice dynamical calculations in the quasiharmonic approximation (QHA) individually for the 
various low-temperature phases (e. g. I, Ha and III in Fig. ^) of LiYF 4 . 

With the advent of high-energy synchrotron X-rays and sophisticated pressure cells, 
high-pressure lattice dynamics has gained considerable momentum in recent years as it 
helps enable a direct comparison of experimental findings, obtained either by the inelastic 
scattering processes or through some modern techniques of ultrafast X-ray diffraction. [11] 
Our calculated variation of phonon frequencies of LiYF 4 in its scheelite phase as a function 
of pressure and temperature (based on a transferable rigid ion model yL M) has already been 
found to auger well with the experimental Raman and infrared data, (jj, |2j 

It is well known that quasiharmonic lattice dynamics (QLD) has the ability to calculate 



12 1 vibrational free energies and other derived properties (e. g. entropy and specific heat) 
of a complex configuration even at very high pressures. To achieve our goal, we make use of 
the same interatomic potential as developed earlier J3( for lattice dynamical calculations. In 
an effort to strengthen the comparison of the three phases of LiYF 4 , phonon dispersions and 
phonon density of states have been examined. Gibbs free energies are then calculated for the 
three high-pressure phases LiYF 4 that are stable at T = 0, with an aim to understand their 
relative stabilities as functions of pressure and temperature in terms of volume compression 
and vibrational entropy. The present work aims at providing a possible pressure-temperature 
(P-T) phase diagram for LiYF 4 to help understand the underlying physics of phase transition 



in such compounds at the microscopic level. The second-order phase transitions and the 
kinetic effects in the first-order phase transitions are specifically investigated by molecular 
dynamics simulations. 

II. VIBRATIONAL PROPERTIES 

A. Raman and infrared active modes 

While the scheelite phase(/4i /a, Z=4) of LiYF 4 has a total of 36 vibrational degrees of 
freedom per primitive cell J2J, La, |5( , the other two high-pressure phases (Ha and III in Fig. |TJ) 
possess 72 phonon modes at each wavevector due to doubling of the primitive cell. Since 
both these high-pressure phases of LiYF 4 belong to the same space group (P2i/c, Z=4), the 
common irreducible representation of the phonons at the zone center is given by 

T: l8A g +l8A u +18B g +18B u , 

of which A g and B 9 modes are Raman active while modes of ungerade symmetry (viz. 
A u and B M ) remain infrared active. It may be noted that the degeneracy of the E g and 
E u modes which were present in phase I (scheelite) of LiYF 4 is lifted in the high-pressure 
phases. Table H] yields a comparison of our calculated zone centre normal modes between 
phases Ha and III at respectively 13 and 18 GPa. Phase Ha has a soft A g mode at 18 
cm -1 which hardens considerably to 92 cm" 1 in phase III. Most other modes harden slightly 
in going from phase Ha to phase III. The vibrational frequencies for each of the phases are 
computed by diagonalization of the respective dynamical matrices using the software DISPR 
[l3j , developed at Trombay. 

B. Phonon dispersion relations 

With a view to compare the phonon dispersion in various phases, we resort to the 
common high symmetry direction of A, which is labeled as the c axis for phase I (I4i/a) 
while it is the b axis for the other high-pressure (monoclinic) phases of similar space group 
symmetry (i. e. P2i/c). On undertaking group theoretical calculations, we obtain for the 
representations of all the normal modes in the three phases as follows: 



Phase I -> A: 8A1 + 8A 2 + 10A 3 (A 3 being doubly degenerate) 
Phase Ha and Phase III -> A: 36Ai + 36A 2 

Phonon dispersions of various polymorphs in LiYF 4 are displayed in Figure El along with 

n n 

the available inelastic neutron scattering data. 2] The experimental data 2] for the phase I 

n 

at zero pressure fit excellently with the lattice dynamical calculations. [3J] Due to doubling 
of the unit cell in going from phase I to phase Ha, the Brillouin zone in phase Ha is halved 
and so there is an apparent folding back of the dispersion branches from zone boundary (of 
phase I) to zone center (of phase Ha). Consequently, there are several optic phonons at the 
zone center in phase Ha at low energies. Interestingly, we come across a very low energy 
zone-center optic mode in phase Ha, indicating the onset of a possible high-pressure phase 
transition at a pressure close to 13 GPa. This espouses our previous MD observations. |8J 
The soft-alike phonon branch of phase Ha may be identified as belonging to the longitudinal 
Kg mode [Table |1] as noted earlier. Further, a crossover between an acoustic phonon mode 
and a low frequency optic mode is observed in phase Ha and also in phase III. We also notice 
that two acoustic phonon branches of A 2 symmetry in the phase III are too closely dispersed 
to remain distinguishable. 

C. Phonon density of states 

To have an overview of the range and extent of various phonon modes in all the three 
phases, we resort to integrating over all the phonons with an energy resolution of 0.5 meV at 
each wave vector on a 16 x 16 x 16 mesh within the irreducible Brillouin zone. As Figure El 
suggests, the energy distribution extends by nearly 10 % in switching from phase I to phase 
Ha while the change in the extent is negligible for the phase III transition. In order to check 
whether it is due to simple pressure effect, we compare the respective frequency distribution 
for phase Ha and phase III at the same pressure of 13 GPa. One may notice that though 
the distribution patterns slightly differ, the extent remain more or less the same. If one 
is further interested in looking for individual contributions of the constituent atoms to the 
entire frequency distribution, partial density of states (PDOS) have to be estimated by 
atomic projections of the one-phonon eigenvectors. 

The labeling of the crystallographic axes viz. a, b, and c varies as per convention 
adopted for specific space groups. In view of these differences and of the need to compare 



the different polymorphs, we choose right-handed orthogonal x, y and z axes as follows: z 
along c and x along a in phase I (J4i/a); z along b and x along c in phases Ha and III 
(P2i/c). Our calculated PDOS for various atoms in various phases are portrayed in Figure^ 
Y atoms are found to contribute largely up to 60 meV while Li atoms make it (especially 
with polarization along x and y) in the higher energy side of the density of states (DOS). 
However, F atom contributions are spread over the entire energy range as if to replicate the 
total DOS to a considerable degree. Such variations in atomic contributions are partly due 
to the mass effect. 

III. THERMODYNAMIC PROPERTIES 

A. Heat capacity and Debye Temperature 

Following the DOS calculations derived out of our simulated data for all the three phases 
of LiYF4, we compare the respective heat capacities at constant pressure(Cp) as a function 
of temperature. FigureG^a) demonstrates how Cp changes for different high-pressure phases. 
We observe that although initially constant pressure heat capacity for phase III runs lower, 
it goes up above the room temperature and in the process even surpasses phase Ha. This 
may be because as we go on increasing temperature (T), the difference Cp(T) - Cy(T) 
[=a v (T)BVT] (where Cy(T) is the constant volume heat capacity and B is the isothermal 
bulk modulus) becomes significantly large, as inset of Figure Efa) suggests, due to the 
existence of much higher thermal expansion ay(T) in phase III than that in phase Ha. A 
detailed discussion on av{T) follows in the next subsection. 

It is illustrative to compare the Debye temperature #d(T) for the various phases [Fig- 
ure|Hfb)] as derived from the calculated Cy(T). We observe that in the very low temperature 
region (say, < 20 K), Op differs a lot from phase to phase (about 14 % increase for I— >IIa 
and about 26 % increase for Ila — ^III) while at higher temperatures (say, 300 K), it comes 
closer for phases Ila and III as difference gets minimal (about 2 % only). 



B. Griineisen parameter and thermal expansion 

Lattice anharmonicity, which leads to a volume dependence of phonon frequencies^), 

is described by the mode Griineisen parameter [lj| given by 

din Ui 
74 = -dln-V (1) 

It is also the only kind of anharmonicity that can be taken care of within the framework of 

quasiharmonic approximation. Moreover, Griineisen parameter is an important quantity as 

it describes the thermoelastic behaviour of materials at high pressures and temperatures. It 

has both the microscopic and macroscopic definition, the former relating to the vibrational 

frequencies of atoms in a material (Eq. QJ) while the latter, to familiar thermodynamic 

properties such as heat capacity and thermal expansion. Unfortunately, the experimental 

determination of 7;, defined in either way, is often not so easy, since the microscopic definition 

requires a detailed knowledge of the phonon dispersion spectrum of a material, whereas the 

macroscopic definition calls for experimental measurements of thermodynamic properties at 

high pressures and temperatures. In this perspective, theoretical model calculation may be 

of great relief. 

Figure El displays our calculated average 7, for various phonon energies and in various 
phases. It can be seen that below 10 meV, there is a significant number of modes with 
negative Griineisen parameter in phase I as well as phase Ha, while phase III modes have 
all positive 7$. 

The effect of pressure on the volume coefficient of thermal expansion(av) can be studied 
through mode Griineisen parameters (7^) in the entire Brillouin zone. In the quasiharmonic 
approximation, each phonon mode of energy Ej contributes to the thermal expansion by way 
of (^y) r YiCvi, where Cyi denotes the specific heat (at constant volume) contribution of the 
jth moc [ e anc [ y the lattice volume. This procedure suits well to the ambient ffuoroscheelite 
system because explicit anharmonicity, which arises mainly out of thermal amplitudes is not 
very significant p| compared to the implicit effect that involves phonon frequency change 
with volume (as obtained from 7^). 

Figure Ufa) demonstrates the variation of ay as a function of temperature. It is also 
apparent that the scheelite as well as the initial high-pressure phase of LiYF 4 have negative 
thermal expansion in the low temperature limit (below 100 K). The reason may be that 



7i has large negative values for phonons below 10 meV in these two phases. However, the 
third phase has no such anomaly, again because it has all positive gammas for phonons of all 
energies (Fig. Ej). Further, the contributions to ay from phonons of different energies (cor- 
responding to various phases of LiYF 4 ) are displayed in Figs. [7(b), (c) and (d) respectively. 
We observe that at low temperatures (e.g. 20 K), contribution of modes upto 10 meV are 
quite significant to ay, but as temperature is increased (e. g. 300 K) higher energy modes 
get more populated and hence, contribute in a large way to the volume thermal expansion. 

C. Mean squared displacements and thermal anisotropy 

In order to gain some insight over how the phonons of various energies are polarized in 
various phases of LiYF4, we plot in Figure|S]tlie partial contributions of these phonons to the 
mean squared thermal amplitude of the constituent atoms. The mean squared displacement 
of atom k along a direction is given by 

oo 

U Qa (k, T) = (u 2 ka ) T = A — [ {n + 1/2) g ka ^)du;, (2) 

7 m k J u 

o 

where n = [exp (^) — l] ; gka{u) is the partial density of states associated with an atom 
k whose mass is m^; A is the normalization constant such that J g ka (u)du> = 1. As FigurelHl 
suggests, the modes at very low energies involve equal displacements of all the atoms that 
correspond to the acoustic modes. Interestingly, between 2 and 9 meV, Y and F atoms have 
larger amplitudes than what relatively light-weight Li atoms possess. It may be noted that 
the basic structure (scheelite) of LiYF4 comprises a pair of strongly bonded LiF4 tetrahedra 
and loosely bonded YF 8 polyhedra 10]. Larger amplitudes of F atoms in the first two phases 
mark the presence of librations of the LiF 4 tetrahedra. 

Our calculated values of equivalent isotropic thermal parameters for Li, Y and F atoms 
in the ambient (i. e. P=0) scheelite phase of LiYF4 are found to be 17, 9 and 14 x 10 _2 A 2 



as against the respective experimental 15] values of 20, 10 and 17 x 10~ 2 A 2 . A detailed 
comparison of the anisotropic thermal parameters (\J aa ) at different temperatures among 
the different phases of LiYF4 is further given in Table |HJ It may be interesting to note that 
while Y and F atoms have comparable atomic displacements along the three directons in 
all the high-pressure phases of LiYF 4 , Li atoms, in contrast, show larger anisotropy along 



the z direction in phase I, along the y direction in phase II and again along the z direction 
in phase III. 

D. Free energy and phase stability 

A P-T phase diagram generally portrays various equilibrium phases at constant tem- 
perature (T) and pressure (P) with the lowest Gibbs free energy (G). From phonon calcu- 
lations, the temperature dependent vibrational free energy at various hydrostatic pressures 
for various phases of LiYF 4 can be estimated in the quasiharmonic approximation. Ther- 
modynamically, we may write [161] 

G n (P,T) = U n -TS n + PV n (3) 

where U n , S n and V n refer respectively to the internal energy, the vibrational entropy and 
the lattice volume of the n th phase. This free energy can be expressed as the sum total 
of configurational contributions (Q c ° n f 19 ) and vibrational contributions(G™ 6 ). The former 
account for mean atomic positions while the latter account for vibrations about the mean 
position. Sometimes, depending on the material to be studied, contributions from magnetic 
effects(G™ a9 ) and electronic effects(G^) can also be significant. However, since LIYF4 is an 
insulator in all the phases being considered, we may use the Born-Oppenheimer approxima- 
tion and therefore we may not explicitly consider the electronic excitations. To repeat it once 
again, the lattice excitations are treated in this work within the quasiharmonic approxima- 
tion where the full Hamiltonian at a given volume is approximated by a harmonic expansion 
about the equilibrium atomic positions, though anharmonic effects are also included to a 
good extent through the implicit volume dependence of the vibrational frequencies. Given 
this, the vibrational Gibbs free energy (G™ fe ) is found to be satisfyingly accurate. 

To include vibrational effects in the present phase diagram, we have calculated dynam- 
ical matrices separately for each of the three phases (viz. I, Ha and III) at pressure intervals 
of 2 GPa on a 4 x 4 x 4 mesh in the irreducible Brillouin zone comprising 64 wave vectors. 
In order to rationalize the behaviour of this dynamical simulation, we calculate 17] the 
enthalpy vs. pressure curves for all the three known structures of LIYF4 and notice that the 
enthalpy changes (due to internal energy and volume) are predominant in pressure- driven 
transitions over free energy changes (due to vibrational energy and entropy). Perhaps this 

8 



is why most of the first-principles phase diagram calculations of pressure-driven phase tran- 
sitions [l8| do not include vibrational entropies, though configurational contributions are 
always taken into account. The other reason for this may be that vibrational entropy dif- 
ferences between phases are assumed to be quite small as we come across in this particular 
case too. 

The stability of a crystalline phase is largely determined by the minimization of the 
Gibbs free energy [l9( and it would, hence, be quite interesting if we can compare the phase- 
wise free energy. The outcome is shown in Figure 01 We note from Fig. Efa) that at 300 
K the free energy plot of the phase I joins smoothly at 6 GPa to that of phase Ha, which 
is consistent with the nature of second order phase transition. Beyond 8 GPa, phase III 
has a lower free energy indicating the greater stability of this phase and also the onset of 
a first order phase transition from phase Ha to phase III. However, the transition pressure 
as obtained through the MD simulations |8|] is higher than 8 GPa due to hysteresis. We 
observe the greater stability of phase III at high pressures arises primarily due to its lower 
volume, while the vibrational energy remains also lower and the entropy becomes higher in 
phase III providing an additional stability. On the plot of differences in vibrational Gibbs 
free energy (AG) as a function of pressure at 300 K, it is clearly seen that G for phase III 
is lower (above 8 GPa) than that for phase Ha. It is also observed (Fig. Efc)) that AG 
between phase Ha and phase III structures decreases with temperature. Another interesting 
point to note from the plots of relative free energy vs. temperature (at 8 GPa) is that in the 
very low temperature region (say, below 100 K), free energy difference shows an anomalous 
behaviour. This can be attributed to the fact that there is a greater density of low frequency 
modes in phase Ha than in phase III due to the soft phonon modes in the former phase. We 
get further support in this regard from Table HI where several vibrational modes have been 
compared between phase Ha and phase III. At low temperatures, these low frequencies are 
populated, thereby lowering the free energy. 

The above comparison of the Gibbs free energies in various phases easily provides the 
phase diagram involving the first order phase transitions. However, the second order tran- 
sitions and also kinetic effects such as hysteresis in the first order transitions are better 
illustrated in molecular dynamics simulations. A comprehensive impression of the phase 
transitions in LIYF4 is given in Fig. El with increasing pressure and temperature (see cap- 



tion of Fig.Hnj). It strengthens our earlier observation J8j that phase II possesses two different 
space groups (viz. phase Ha: P2\jc ; phase lib: 12/ a) existing in two different ranges of 
temperaure. The low temperature phase Ha involves small displacements of some of the F 
atoms from their ideal positions of body-centered symmetry in phase lib. The transition 
Ha — ► III occurs at a much higher pressure of 14 GPa than the equilibrium phase boundary 
at 8 GPa due to hysteresis at the time scale of the simulation and is found irreversible on 
release of pressure. 

IV. CONCLUSIONS 

We have been able to demonstrate how a conjugate lattice dynamics - molecular 
dynamics study can eventually lead one to calculate the vibrational free energy and other 
essential thermodynamic functions of LiYF 4 in all its available high-pressure phases from 
the mere knowledge of an interatomic pair potential. It also sheds light on the useful under- 
standing of some important physical properties and phenomena associated with the phase 
equilibria. Finally, our calculated results have the potential ingredients to help analyze the 
experimental inelastic scattering data, if high-pressure (along with low temperature) phonon 
measurements are carried out in future, which could be affirmative as well as interesting too. 
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TABLE I: Comparison of calculated zone-center phonon frequencies between phase Ha (P= 13 
GPa) and phase III (P=18 GPa) of LiYF 4 at T=0 K. 
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TABLE II: Calculated anisotropic thermal parameters (in units of 1CP 4 A 2 ) for the constituent 
atoms at various temperatures of 20 and 300 K associated with the three high-pressure phases of 
LiYF,4. It may be noted that in the scheelite phase (J4i/a, Z=4), all the F atoms are symmetrically 
equivalent. See text for the labeling of x, y, z directions. 
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Figure Captions 

FIG. 1: Crystal structures of LiYF4 belonging to (a) phase I, (b) phase Ha, (c) phase lib and (d) 
phase III. 

FIG. 2: Comparison of phonon dispersion relations along the same high symmetry direction of A 
in the three phases of LIYF4 at respectively 0, 13 and 18 GPa. Lines refer to our calculations while 
symbols represent the experimental data 



a. 



FIG. 3: (a)Calculated phonon density of states, g(E), in the three phases of LIYF4 at respectively 
4, 13 and 18 GPa; (b) comparison of g(E) in the two high-pressure phases at the same pressure of 
13 GPa. 

FIG. 4: Calculated partial density of states for the constituent atoms with contributions from 
polarizations along three orthogonal directions in the three phases of LIYF4 at respectively 4, 13 
and 18 GPa. Since phase I is of tetragonal symmetry, the contributions along x and y directions are 
identical for this phase. The relationship between the orthogonal directions and the crystallographic 
axes is given in section IIC. 

FIG. 5: Calculated (a) heat capacity (Cp) and (b) Debye temperature^/?) as a function of tem- 
perature for the three phases of LIYF4 at respectively 4, 13 and 18 GPa. The inset depicts the 
difference of Cp(T) and Cy(T) over the same temperature region. 

FIG. 6: Calculated average mode Griineisen parameter (ji) as a function of phonon energy for the 
three phases of LIYF4 at respectively 4, 13 and 18 GPa. 

FIG. 7: Calculated volume thermal expansion coefficient (ay) as a function of temperature for the 
three phases of LIYF4 at respectively 4, 13 and 18 GPa. (b), (c) and (d) demonstrate how phonons 
of different energies (in phase I, phase Ha and phase III respectively) contribute to ay at 20 and 
300 K. 

FIG. 8: Calculated thermal mean squared displacements (on a semi-logarthmic scale) of the con- 
stituent atoms for (a) phase I, (b) phase Ha and (c) phase III due to phonons of various energies 
at 300 K. Insets in (a), (b) and (c) refer to the same but on an expanded scale up to 10 meV for 
the respective phases of LIYF4 at 4, 13 and 18 GPa. It may be noted that in phases Ha and III, 
there are four symmetrically different F atoms, however, only the average contributions of all the 
F atoms are shown here. 
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